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Abstract 

, '•y ' Many applications in network analysis require algorithms to sample uniformly at random from 

I the set of all graphs with a prescribed degree sequence. We present a Markov chain based approach 

which converges to the uniform distribution of all realizations for both the directed and undirected 
case. It remains an open challenge whether these Markov chains are rapidly mixing. 
^ For the case of directed graphs, we also explain in this paper that a popular switching algorithm 

■ fails in general to sample uniformly at random because the state graph of the Markov chain decom- 

' poses into different isomorphic components. We call degree sequences for which the state graph is 

, strongly connected arc swap sequences. To handle arbitrary degree sequences, we develop two dif- 

^ ' ferent solutions. The first uses an additional operation (a reorientation of induced directed 3-cycles) 

^ ' which makes the state graph strongly connected, the second selects randomly one of the isomor- 

1 I , phic components and samples inside it. Our main contribution is a precise characterization of arc 

swap sequences, leading to an efficient recognition algorithm. Finally, we point out some interesting 
I consequences for network analysis. 

in : 

^ ; 1 Introduction 

• ' We consider the problem of sampling uniformly at random from the set of all realizations of a prescribed 

I . degree sequence as simple, labeled graphs or digraphs, respectively, without loops. 

o^: 

' Motivation. In complex network analysis, one is interested in studying certain network properties of 

^ I some observed real graph in comparison with an ensemble of graphs with the same degree sequence 

r'"! ■ to detect deviations from randomness |MKI+04| . For example, this is used to study the motif content 

rS I of classes of networks lMIK+041 . To perform such an analysis, a uniform sampling from the set of all 

' realizations is required. A general method to sample random elements from some set of objects is via 

rapidly mixing Markov chains |Sin92i ISin93| . Every Markov chain can be viewed as a random walk on a 
directed graph, the so-called state graph. In our context, its vertices (the states) correspond one-to-one 
to the set of all realizations of prescribed degree sequences. For a survey on random walks, we refer to 
Lovasz |Lov96| . 

A popular variant of the Markov chain approach to sample among such realizations is the so-called 
switching- algorithm. It starts with a given realization, and then performs a sequence of 2-swaps. 

In the undirected case, a 2-swap replaces two non-adjacent edges {a, 6}, {c, d} either by {a, c}, {6, d} 
or by {a, d}, {6, c}, provided that both new edges have not been contained in the graph before the 
swap operation. Likewise, in the directed case, given two arcs (a, &), (c, d) with all vertices a, 6, c, d being 
distinct, a 2-swap replaces these two arcs by (a, d), (c, b) which are currently not included in the realization 
(the latter is crucial to avoid parallel arcs). The switching algorithm is usually stopped heuristically after 
a certain number of iterations, and then outputs the resulting realization as a "random element". 



*This work was partially supported by the DFG Focus Program Algorithm Engineering, grant Mu 1482/4-1, and by a 
VolkswagenStiftung grant for the project "Impact on motif content on dynamic function of complex networks". 
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Figure 1: Example of digraphs where no 2-swap operation can be applied. 

For undirected graphs, one can prove that this switching algorithm converges to a random stage. The 
directed case, however, turns out to be much more difficult. The following example demonstrates that 
the switching algorithm does not even converge to a random stage. 

Example 1.1. Consider the following class of digraphs D ~ (V, A) with3n verticesV ~ {wi, W2, . • ■ , W3„}, 
see Figure{^ Roughly speaking, this class consists of induced directed 3-cycles Ci formed by triples Vt = 
{v3i,V3i+i,V3i+2} of vertices, and arcs Ai = {(113^, wsi+i), (113^+1, W3i+2), (w3i+2, t'sOl « G {0,...,n-l}. 
All vertices of cycle C'l are connected to all other vertices of cycles with larger index than i. More formally, 
let A' := {{v,w)\v G V^,w G Vj,i < j}. We set A := A' U (Uf^i A;). 

It is easy to check that no 2-swap can be applied to this digraph. However, we can independently 
reorient each of the n induced 3-cycles, leading to 2"/'^ many (isomorphic) realizations of the same degree 
sequence. Thus, if we use a random walk on the state graph of all realizations of this degree sequence and 
use only 2-swaps to define the possible transitions between realizations, this state graph consists exactly 
of 2"/^ many singleton components. Hence, a "random walk" on this graph will be stuck in a single 
realization although exponentially many realizations exist. 

More examples of graph classes of this type will be given in the Appendix. It is interesting to note 
that 2-swap operations suffice to sample directed graphs with loops as has been proven by Ryser |Rys57| 
in the context of square matrices with {0, l}-entries which can be interpreted as node-node adjacency 
matrices of digraphs with loops. 

Realizability of degree sequences. In order to use a Markov chain approach one needs at least 
one feasible realization. In applications from complex network analysis, one can usually take the degree 
sequence of some observed real world graph. Otherwise, one has to construct a realization. 

The realization problem, i.e., characterizing the existence and finding at least one realization, has 
quite a long history. First results go back to the seminal work by Tutte who solved the more general 
/-factor problem |Tut52| . Given a simple graph G = {V,E) and a function / : V{G) 1—)- No, an f -factor 
is a subgraph H oi G such that every vertex z; G ^ in this subgraph H has exactly degree dciv) = f{v). 
Tutte gave a polynomial time transformation of the /-factor problem to the perfect matching problem. 
This implies the first polynomial time algorithm for finding some /-factor |Tut54| . For a survey on 
efficient algorithms for the /-factor problem by matching or network flow techniques, we refer to Chapter 
21 of Schrijver |Sch03| . Clearly, if the given graph G is complete, then every /-factor is a solution of 
the degree sequence problem. Erdos and Gallai |EG60| proved a simpler Tutte-type result for the degree 
sequence problem. Already in 1955, Havel |Hav55| developed a simple greedy-like algorithm to construct 
a realization of a given degree sequence as a simple undirected graph without loops. A few years later, 
Hakimi |Hak62l IHak65| studied the simpler case of undirected graphs with multiple edges. 

It is also well-known how to test whether a prescribed degree sequence can be realized as a digraph. 
Chen |Che66| presented necessary and sufficient conditions for the realizability of degree sequences which 
can be checked in linear time. Again, the construction of a concrete realization is equivalent to an /-factor 
problem on a corresponding undirected bipartite graph. Kleitman and Wang |KW73| found a greedy- type 
algorithm generalizing previous work by Havel [Hav55] and Hakimi |Hak62i IHak65| . This approach has 
recently been rediscovered by Erdos et al. [EMT09) . 

Related work. Kannan et al. |KTV99| showed how to sample bipartite undirected graphs via Markov 
chains. They proved polynomial mixing time for regular and near-regular graphs. Cooper et al. |CDG07] 
extended this work to non-bipartite undirected, d-regular graphs and proved a polynomial mixing time 
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for the switching algorithm. More precisely, they upper bounded the mixing time in these cases by 
d}^n^{dn\og{dn) + log(£^^)), for graphs with \V\ = n. In a break-through paper, Jerrum, Sinclair, and 
Vigoda |JSV04| presented a polynomial-time almost uniform sampling algorithm for perfect matchings in 
bipartite graphs. Their approach can be used to sample arbitrary bipartite graphs and arbitrary digraphs 
with a specified degree sequence in 0(n^^ log^ n) via the above-mentioned reduction due to Tutte. In the 
context of sampling binary contingency tables, Bezakova et al. |BBV07| managed to improve the running 
time for these sampling problems to 0{n^^ log^ n), which is still far from practical. 

McKay and Wormald |MW90[ IMW91| use a configuration model and generate random undirected 
graphs with degrees bounded by o(n^/^) with uniform distribution in 0{m^dmax) time, where dmax 
denotes the maximum degree, and m the number of edges. Steger and Wormald |SW99| introduced a 
modification of the configuration model that leads to a fast algorithm and samples asymptotically uniform 
for degrees up to o{n^^^^). Kim and Vu |KV03| improved the analysis of Steger and Wormald's algorithm, 
proving that the output is asymptotically uniform for degrees up to 0{n^^^~^), for any e > 0. Bayati 
et al. [BKS09| recently presented a nearly-linear time algorithm for counting and randomly generating 
almost uniformly simple undirected graphs with a given degree sequence where the maximum degree is 
restricted to d^ax — 0(m^/'*~'^), and r is any positive constant. 

Random walks and Markov chains. Let us briefly review the basic notions of random walks and 
their relation to Markov chains. See |Lov96[ IJS96I ISin93| for more details. A random walk (Markov 
chain) on a digraph D = (V, A) is a sequence of vertices wq, wi, . . . , ut, . . . where (ui, Wi+i) G A. Vertex Vq 
represents the initial state. Denote by the out-degree of vertex v gV. At the tth step we move to 

an arbitrary neighbor oi Vt with probability l/d~jj{vt) or stay at Vt with probability (1 — 1/(111)) /d~jj{vt), 
where I'ivt) denotes the number of neighbors of vt- Furthermore, we define the distribution of V at time 
t G Z+ as the function Pt G [0, with Pt{i) := Prob{vt = i). A well-known result |Lov96| is that Pt 
tends to the uniform stationary distribution for f — !> cx), if the digraph is (1) non-bipartite (that means 
aperiodic), (2) strongly connected (i.e., irreducible), (3) symmetric, and (4) regular. A digraph D is 
do -regular if all vertices have the same in- and out-degrees do- 

In this paper, we will view all Markov chains as random walks on symmetric d^-regular digraphs 
D ~ (V, A) whose vertices correspond to the state space V . The transition probability on each arc 
(u, w) & A will be the constant l/dn. 

Our contribution. In this paper, we prove the following results. 

• For undirected graphs we analyze the well-known switching algorithm. It is straight-forward to 
translate the switching algorithm into a random walk on an appropriately defined Markov chain. 
This Markov chain corresponds to a symmetric, regular, strongly connected, non-bipartite simple 
digraph with directed loops allowed. Thus, it converges to the uniform distribution of all realiza- 
tions. Each realization of the degree sequence is a vertex of this digraph, and two realizations are 
mutually connected by arcs if and only if their symmetric difference is an alternating 4-cycle (i.e., 
corresponds to a 2-swap). This graph becomes regular by adding additional loops, see Section [51 

Cooper et al. |CDG07| already showed in the context of regular graphs that the underlying digraph 
of this Markov chain is strongly connected, but we give a much simpler proof of this property. Its 
diameter is bounded by the number m of edges in the prescribed degree sequence. 

• Carefully looking at our Example II. IL we observe that in the directed case the state graph becomes 
strongly connected if we add a second type of operation to transform one realization into another: 
Simply reorient the arcs of an induced directed 3-cycle. We call this operation 3-cycle reorientation. 
We give a graph-theoretical proof that 2-swaps and 3-cycle reorientations suffice not only here, but 
also in general for arbitrary prescribed degree sequences. These observations allow us to define 
a Markov chain, very similar to the undirected case. The difference is that two realizations are 
mutually connected by arcs if and only if their symmetric difference is either an alternating directed 
4-cycle or 6-cycle with exactly three different vertices. Again, this digraph becomes regular by 
adding additional loops, see Section [31 The transition probabilities are of order 0{l/m^), and the 
diameter can be bounded by 0(rn), where m denotes the number of arcs in the prescribed degree 
sequence. 
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In the context of (0, l)-matrices with given marginals (i.e., prescribed degree sequences in our 
terminology) , Rao et al. |RJB96| similarly observed that switching operations on so-called "compact 
alternating hexagons" are necessary. A compact alternating hexagon is a 3 x 3-submatrix, which 
can be interpreted as the adjacency matrix of a directed 3-cycle subgraph. They define a random 
walk on a series of digraphs, starting with a non-regular state graph which is iteratively updated 
towards regularity, i.e. their Markov chain converges asymptotically to the uniform distribution. 
However, it is unclear how fast this process converges and whether this is more efficient than starting 
directly with a single regular state graph. Since Rao et al. work directly on matrices, their transition 
probabilities are of order 0{\/n^), i.e., by several orders smaller than in our version. 

Very recently, Erdos et al. |EMT09| proposed a similar Markov chain approach using 2-swaps and 
3-swaps. The latter type of operation exchanges a simple directed 3-path or 3-cycle (ui, V2), (w2, ^3), 
(^3: (the first and last vertex may be identical) by (vi, U3), (W3, V2), (w2, V4), but is a much larger 
set of operations than ours. 

• Although in directed graphs 2-swaps alone do not suffice to sample uniformly in general, the cor- 
responding approach is still frequently used in network analysis. One reason for the popularity 
of this approach — in addition to its simplicity — might be that it empirically worked in many 
cases quite well |MKI+04| . In this paper, we study under which conditions this approach can be 
applied and provably leads to correct uniform sampling. We call such degree sequences arc-swap 
sequences, and give a graph-theoretical characterization which can be checked in polynomial time. 
More specifically, we can recognize arc-swap sequences in O(m^) time using matching techniques. 
Using a parallel Havel-Hakimi algorithm by LaMar |LaM09| . originally developed to realize Euler 
sequences with an odd number of arcs, the recognition problem can even be solved in linear time. 
This algorithm also allows us to determine the number of induced directed 3-cycles which appear 
in every realization. 

However, the simpler approach comes with a price: our bound on the diameter of the state graph 
becomes mn and so is by one order of n worse in comparison with using 2-swaps and 3-cycle 
reorientations. Since half of the diameter is a trivial lower bound on the mixing time and the 
diameter also appears as a factor in known upper bounds, we conjecture that the classical switching 
algorithm requires a mixing time with an order of n more steps as the variant with 3-cycle 
reorientation. 

In those cases where 2-swaps do not suffice to sample uniformly, the state graph decomposes into 2^ 
strongly connected components, where k is the number of induced directed 3-cycles which appear in 
every realization. We can also efficiently determine the number of strongly connected components of 
the state graph (of course, without explicitly constructing this exponentially sized graph). However, 
all these components are isomorphic. This can be exploited as follows: For a non-arc-swap sequence, 
we first determine all those induced directed 3-cycles which appear in every realization. By reducing 
the in- and out-degrees for all vertices of these 3-cycles by one, we then obtain a new sequence, now 
guaranteed to be an arc-swap sequence. On the latter we can either use the switching algorithm or 
our variant with additional 3-cycle reorientations on a smaller state graph with a reduced diameter 
n(m — 3fc) or m — 3fc, respectively, yielding an important practical advantage. 

Our results give a theoretical foundation to compute certain network characteristics on unlabeled 
digraphs in a single component using 2-swaps only. For example, this includes the analysis of 
the motif content |MSOI+02| . Likewise we can still compute the average diameter among all 
realizations if we work in a single component. However, for other network characteristics, for 
example betweenness centrality on edges |KLP"'"05| . this leads in general to incorrect estimations. 

Overview. The remainder of the paper is structured as follows. In Section [5] we start with the undi- 
rected case. We introduce appropriately defined state graphs underlying our Markov chains, and show 
for these graphs crucial properties like regularity and strong connectivity. We also upper bound their 
diameter. The more difficult directed case is presented in Section [31 Afterwards, in Section 2) we char- 
acterize those degree sequences for which a simpler Markov chain based on 2-swaps provably leads to 
uniform sampling in the directed case. We also describe a few consequences and applications. Finally, 
we conclude with a short summary and remarks on future work. 
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2 Sampling Undirected Graphs 



In this section we show how to sample undirected graphs with a prescribed degree sequence uniformly 
at random with a random walk. This section is structured as follows. We first give a formal problem 
definition and introduce some notation. Then we introduce an appropriately defined Markov chain and 
prove that it has all desired properties. 

Formal problem definition. In the undirected case, a degree sequence S of order n is the ordered 
set (ai, 02, . . . , a„) with a; G Z+,ai > 0. Let G = {V,E) be an undirected labeled graph G = {V,E) 
without loops and parallel edges and \V\ ~ n. We define the degree-function d : V ^ Z+ which assigns 
to each vertex Vi € V the number of incident edges. We call S a graphical sequence if and only if there 
exists at least one imdirected labeled graph G = (V, E) without any loops or parallel edges which satisfies 
d{vi) = Gi for all Ui G V and i G {1, . . . , Any such undirected graph G is called realization of S. 

We define an alternating walk P for a graph G ~ {V, E) as a sequence P := (vi ,V2, ■ ■ ■ ,vg) of vertices 
Vi & V where either {vi,Vi+i} G E{G) and {dj_i,Wi} ^ E{G) or {vi,Vi+i} ^ E{G) and {vi-i,Vi] G E{G) 
for i mod 2=1. The length of a walk (or path, cycle, respectively) is the number of its edges. We call an 
alternating walk C of even length alternating cycle if vi = vg is fulfilled. For two realizations G,G', the 
symmetric difference of their edge sets E{G) and E{G') is denoted as GAG" := {E{G)\E{G')) U {E{G') \ 
E{G)). A graph is called Eulerian if every vertex has even degree. Note that the symmetric difference 
GAG' of two realizations G, G' is Eulerian and hence always decomposes into a number of alternating 
cycles. 

The Markov ciiain. We denote hy ^ = {V[f,,A^) the digraph for our random walk, the state graph, 
for short. Its underlying vertex set Vjf, is the set of all realizations of a given degree sequence S. For a 
realization G, we denote by Vq the corresponding vertex in VJ/,. The arc set is defined as follows. 

a) We connect two vertices Vg,Vg' G V^,G ^ G' with arcs (Vg,Vg') and (Vg',Vg) if and only if 
I GAG" I = 4 is fulfilled. 

b) We set for each pair of non-adjacent edges {vi-^^Vi^}, {vi^,Vi^} G E{G),ij G {l,...,n} a directed 
loop iVG,VG) if and only if {v,„v,,} G E{G) V {v,„v^,} G E{G). 

c) We set for each pair of non-adjacent edges {wji , Wia}, {wi.,, Wii} G E{G),ij G {l,...,n} a directed 
loop {Vg,Vg) if and only if {v,,,v,,} G E{G) V v.J G E{G). 

d) We set one directed loop (Vg, Vg) for all Vg &V^. 

Lemma 2.1. T/ie state graph ^ = (V^^, A^,) is non-bipartite, symmetric, and regular. 

Proof. Non-bipartiteness follows from the insertion of directed loops. Likewise, symmetry is obvious since 
we always introduce arcs in both directions in case a). For each pair of non-adjacent edges of a realization 
G, we introduce exactly two arcs in ^. These arcs either connect two neighboring states or are directed 
loops. Thus each vertex Vg & V^p has an out-degree of twice the number of non-adjacent edges in G plus 
one (for the loop in step d)). Due to symmetry, the out-degree equals the in-degree. For each realization 

G, the number of pairs of non-adjacent edges is exactly ('^') — J2vi&viG) ('^'^2"'^) ~ Cf ) ~ (2')' 
that is a constant independent of G. □ 

The next step is to show that the state graph is strongly connected. We first prove the following 
auxiliary proposition which asserts that the symmetric difference of two different realizations always 
contains a vertex-disjoint path of length three. 

Proposition 2.2. Let S be a graphical sequence and G and G' be two different realizations, i.e., GAG' ^ 
0. Then there exists a vertex-disjoint alternating walk P = (vi, V2, W3, W4) in GAG' with {vi, V2}, {ws, U4} G 
E{G) and {w2,U3} e E{G'). 

Proof. In the proof of this proposition, we argue only about edges in the symmetric difference GAG' 
which is assumed to be non-empty. Therefore, there are edges {vi, V2}, {v2, V3} with {vi, V2} G E{G) and 
{^'2,^3} G E{G') and vi ^ V3. If there is also an edge {v3,V4} G E{G) with V4 ^ vi, we are done with 
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the vertex-disjohit alternating walk P ~ (fi, ^2, W3, 'y4) as desired. Otherwise, the symmetric difference 
must contain the edge {vz,vi} € E{G) and also some edge {wi,W4} G E{G'). Note that W4 7^ V2 and 
va v-i. This implies the existence of another edge {W4, w.g} g E{G). Note also that 7^ W3, since we are 
in the case that {v^, v^} does not exist. Either U5 = V2 or w.g is a new vertex disjoint from {vi, . . . , W4}. 
Therefore, in both cases P = (173, ui, U4, U5) is a vertex-disjoint alternating walk composed of edges from 
the symmetric difference. □ 

Lemma 2.3. Let S be a graphical sequence and let G ^ G' be two realizations. Then there exist realiza- 
tions G'o,Gi, ...,Gk with Go := G, Gk :== G' and \GiAG^+i\ = 4 where k < ^\GAG'\ - 1. 

Proof. We prove the lemma by induction according to the cardinality of the symmetric difference |GAG' | = 
2k. For K := 2 we get |GAG'| = 4. The correctness of our claim follows with Gi := G'. We assume the 
correctness of our claim for all k < i. Consider |GAG'| — 2i + 2. According to Proposition 12.21 there 
exists in GAG' an alternating vertex-disjoint walk P = (wi, W2, "^3, Vi) with {vi, V2}, {^3, V4} £ E{G) and 
{V3,V2} e E{G'). 

case 1: Assume {vi,Vi} g E{G') \ E{G). 

This implies {vi, v^} g GAG'. Gi := (Go \ {{vi, V2}, {v^, v^}}) U {{^2, ^'3}, {vi, ^4}} is a realization 
of S and it follows |GoAGi| = 4 and |GiAG'| = 2£ + 2 - 4 = 2(^ - 1). Note that after this 
step, GiAG' may consist of several connected components, but each of them has strictly smaller 
cardinality. Thus, we obtain realizations Gi, G2, . . . ,Gk with Gk ■= G' and |GiAGi+i| = 4 where 
k-1 < i|GiAG'| - 1. Hence, we get the sequence Go,Gi, . . . , Gfc with fc = 1 + i|GiAG'| - 1 = 
i(|GAG'| -4) < i|GAG'| - 1. 

case 2: Assume {vi,Vi} g E{G) r\E{G'). 

This implies {vi, v^} ^ GAG' . P is an alternating subpath of an alternating cycle G ~ (u4, Vi, . . . , vj, 
vi, V2, W3, W4) with {vi, W4}, {vi, Vj} g A(G') of length |G| > 6. We construct a new alternating cycle 
G* :— {G\P)U{{vi, V4}} with length |G* | = |G|-2 > 4. We swap the arcs in G* and get a reahzation 
G* of S with |GoAG*| = |G*| < 2€ and |G*AG'| = |GAG'| - (|G*| - 1) -h 1 < 2^. The symmetric 
difference G*AG' may consist of several connected components, but their total length is bounded 
by 2£. Thus there exist sequences Gj := G, Gj, . . . , G^^ := G* and Gg := G*, G?, . . . , Gl^ = G' 
with fci < i|GoAG*|-l = i|G*|-l and k2 < i|G*AG'|-l = i(|GAG'| - (|C*| - 1) + 1) - 1 < 21 
We arrange these sequences one after another and get a sequence which fulfills k = ki + k2 = 
i|C*| - 1 + i(|GAG'| - (|G*| - 1) + 1) - 1 = i(|GAG'|) - 1. 

case 3: Assume {vi,Vi} g E{G) \E{G'). 

This implies {vi,V4} g GAG'. Assume first that the symmetric difference GAG' contains an al- 
ternating cycle G which avoids P. Then, we can apply the induction hypothesis to G. Swap- 
ping the edges of G, we get a realization G* of sequence S with |GoAG*| = |G*| < 2i and 
G*AG'| = I GAG' I — |G| < 2£. According to the induction hypothesis there exist sequences 
Gl := G,G\,...,Gl^ G* and G^ := G*, G^, . . . , G^^ = G' with fci < i|GoAG*| - 1 and 
^2 < ^|G*AG'| — 1. We arrange these sequences one after another and get a sequence which fulfills 
fc = fci + fc2 < 5|GoAG*| - 1 + i|G*AG'| - 1 < i(|GAG'|) - 1. 

It remains to consider the case that such a cycle C does not exist. In other words, every alternating 
cycle in GAG' includes edges from P. 

The alternating walk P can be extended to an alternating cycle C* ~ (wi, U2, W3, ^'4, ws, • . • , W2t, ^i), 
t > 3 using only arcs from GAG'. To construct G*, start with P, and keep adding alternating 
edges until you reach the start vertex vi for the first time with an edge {vi,vi} g E{G'). Since 
the symmetric difference is Eulerian, you will not get stuck before reaching vi with such an edge. 
Note that G* must contain the edge {wi,W4}, as otherwise an alternating cycle of type G would 
exist. This also implies the existence of an alternating sub-walk Pi = (114, W5, . . . ^vq^v^) of G* of 
odd length (at least of length 3), starting and ending with edges in E{C*). Likewise, there must 
be another alternating sub-walk P2 = {wi, U7, . . . , us, wi}, also of odd length (at least of length 3), 
starting and ending with edges in E{C*). The situation is visualized in Figure |5] 

In this scenario, we have v^, ^ vj, as otherwise {E{Pi) \ {{vj, wi}}) U {{w? = W5, W4}, {^4, wi}} would 
be an alternating cycle of the form we have excluded above. We have four subcases with respect to 
the existence of edges between ^5 and W7. 
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edges in G 
edges in G' 



Figure 2: Proof of Lemma [2.31 Edges of the symmetric difference GAG' in case 3. 
case a) {vg, V7} G E{G) \ E{G'): 

This would imply the existence of the alternating cycle C = (?J5, U4, ui, W7, W5), excluded above, 
case b) {V5,V7} e E{G')\E{Gy. 

This would imply the existence of the alternating cycle C = (vj, 1)5, . . . ,ve,V4,Vi,Vs, ■ ■ ■ jV^), also 
excluded above. 

case c) {w5, V7} £ E{G) n E{G'): 

Then there is an alternating cycle C" = (^7, V5, U4, 1'l, uy) on which we can swap the edges in a single 
step. This leads to a realization G* = d with IGqAG*! = 4 and |G*AG'| = 2£. 

case d) {wg, W7} ^ -B(G) and {v^^vr} E(G'): 

As in case c), we consider the alternating cycle G' = (wy, W5, Pi\{{t;4, ws}}, {wi, 'y4}, -P2\{{i'7, wi}}, uy). 
Swapping edges on G', we get a realization G* = Gk-i, but this time. |GoAG*| = |G'| < 2^ and 
|G*AG'| ^ IGAG'I - (|G'| - 1) + 1 < 2£. According to the induction hypothesis there exist se- 
quences Gj := G,Gl,...,Gl^ := G* and Gg := G* ,Gl, . . . ,Gl^ ^ G' with fci < i|GoAG*| - 1 = 
i(|GAG'|-|C'l+2)-landA:2 < i|G*AG'|-l = We arrange these sequences one after an- 

other and get a sequence which fulfills = fci+/s2 < i(|GAG'|-|C'l+2)-l+5|G'|-l = \\GAG'\-l. 

case 4: Assume {ui,U4} ^ £^(G) U E{G'). This implies {wi,U4} ^ GAG'. It exists the alternating cycle 
G := {vi,V2,V3,V4,vi)) with {wi,?;4} ^ E{G). Gi := (Go \ {{ui, W2}, {fs, W4}}) U {{^3, 'y2}, {^'i, ""4}} 
is a realization of 5 and it follows |GoAGi| = 4 and |GiAG'| = 2£ + 2 - 2 = 2^. According to 
the induction hypothesis there exist realizations Gi,G2,...,Gfc with Gk ■= G' where and k < 
i|GAG'| - 2. Hence, we get the sequence Go,Gi,...,Gfc with k < i|GAG'| - 1. 

□ 

We have shown that the state graph = {V^,A^) is a d-regular, symmetric, non- bipartite, and 
strongly connected digraph. Hence, the corresponding Markov chain has the uniform distribution as 
its stationary distribution. A random walk on = can be described by Algorithm [1] This 

algorithm requires a data structure DS containing all pairs of non-adjacent edges in G. 



3 Sampling Digraphs 

We now turn the directed case. As before, we start with the formal problem definition and some additional 
notation. Then, we introduce our Markov chain and analyze its properties. 

Formal problem definition In the directed case, we define a degree sequence S" as a sequence of 
2-tuples (^(11), (11),. . . , (12)^ with a„6, e Z+ ,i e {1, . . . ,n} where Oi > or 6, > 0. 

Let G = {V,A) be a directed labeled graph G = {V,A) without loops and parallel arcs and \V\ = n. 
We define the in- degree- function : V ^ which assigns to each vertex Vi £ V the number of 
incoming arcs and the out- degree- function d^'.V^ which assigns to each vertex Vi E V the number 
of outgoing arcs. We denote S as graphical sequence if and only if there exists at least one directed labeled 
graph G = {V,A) without any loops or parallel arcs which satisfies dg(t;j) = bi and d^lvi) = a.; for all 
Vi E V and i G {1, . . . , n}. Any such graph G is called realization of S. Let iJ be a subdigraph of G. We 
say that H = {Vh,Ah) is an induced subdigraph of G if every arc of A with both end vertices in Vh is 
also in Ah- We write H = G {Vh) ■ 
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Algorithm 1 Switching Algorithm 

Input: sequence S, an undirected graph G = (V, E) with dcivi) ~ ai for all i G {1, . . . , n} and Vi G V , 
a mixing time r. 

Output: A sampled undirected graph G" = (y, E') with dcivi) — ai for all i G {1, . . . , n} and G V. 



1: t := 0, G' := G / /initialization 

2: while t < T do 

3: Choose an element p from DS uniformly at random./ /p is a pair of non-adjacent edges. 

4: Let p be the pair of edges {vi^ , }, {vi, ,Vi^}. 

5: Choose with probability i between case a) and case b). 

6: if case a) then 

7: if {vii,ViJ,{vi^,Vi^} ^ E{G') then 

8: // Either walk on to an adjacent realization 

9: Delete {vi^,Vi^},{vi^,Vi^} in E{G'). 

10: Add {w,i,w,4},{wj3,wjj to E{G'). 

11: else 

12: 11 or walk a loop: 'Do nothing' 

13: end if 

14: else 

15: / / case b) 

16: if {vi^ , ■Uig } , {vi., ,Vi^} i E{G') then 

17: / / Either walk on to an adjacent realization 

18: Delete {v^^,Vi^},{vi^,Vi^} in E{G'). 

19: Add {v^^,v^.,},{v,^,v,^} to E{G'). 

20: else 

21: 11 or walk a loop: 'Do nothing' 

22: end if 

23: end if 

24: update data structure DS 

25: t ■i^t + l 

26: end while 



The symmetric difference GAG' of two realizations G ^ G' \& defined analogously to the undirected 
case. Consider for example the realizations G and G' with A{G) := {(wi, W2), (ws, W4)} and A{G') := 
{(ui, U4), (y3, V2)} consisting of exactly two arcs. Then the symmetric difference is the alternating directed 
4-cycle C := {vi,V2,V3,V4,vi) where {vi,Vi+i) G A{G) for i G {1,3} and {vi+i,Vi) G A{G') taking indices 
i mod 4. We define an alternating directed walk P for a directed graph G = (V, A) as a sequence P := 
{vi,V2, ■ ■ ■ ,vi) of vertices Vi G F where either {vi,Vi+i) G A{G) and ^ ^(G) or {vi,Vi+i) ^ A(G) 

and {vi,Vi-i) G ^(G) for i mod 2 = 1. We call an even alternating directed walk G alternating directed 
cycle if vi — vi is fulfilled. The symmetric difference of two realizations always decomposes into a number 
of alternating directed cycles, see Figs. [3] and HI 

The Markov chain. In the directed case, we denote the state graph for our random walk by $ = 
{V^,A^). Its underlying vertex set is the set of all realizations of a given degree sequence S. For a 
realization G, we denote by Vq the corresponding vertex in V^. The arc set A^ is defined as follows. 

a) We connect two vertices Vg, Vq' G V^, G ^ G' with arcs (Vg, Vg') and (Vg", Vq) if and only if one 
of the two following constraints is fulfilled 

1. IGAG'I =4 

2. I GAG' I = 6 and GAG' contains exactly three different vertices. 

b) We set a directed loop (Vg, Vg) 

1. for each pair of non-adjacent arcs {vi-^^,Vi^),{vi^,Vi^) G A{G)^ij G {l,...,n} if and only if 
{vinVi^) G A{G) V (uigjWij) G A{G) in a realization G, 
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^► inG 
♦ inG' 



Figure 3: Example: Two realizations G and G' . 

(2)^--^.@^ (5 

C ^ 




Figure 4: Decomposition of the symmetric difference GAG" of Fig. [3] into a minimum number of alter- 
nating directed cycles. 



2. for each directed 2-path (w^^ , ), (wij , ) G ^(G) if and only if one of the following constraints 
is true for a realization G, 

i) G ^(G) V (w.3,w»J e A(G) V K,w,;3) e A(G), 

iii) 13 < ii V 13 < 12. 

3. if G contains no directed 2-path. 

Lemma 3.1. The state graph $ := (V^, A^) is non-bipartite, symmetric, and regular. 

Proof. In our setting we connect two vertices at each time in both directions. Hence. $ is symmetric. Fur- 
thermore, if some realization G contains no directed 2-path, then each G is a realization of a sequence 5*, 
only consisting of sinks and sources. With our setting $ contains for each Vq G a directed loop and is 
therefore non-bipartite, see item 6)3 in our construction. Let us now assume that a realization G contains 
a directed 2-path. Either there exists a third arc which completes these two arcs to a directed 3-cycle 
or not. In all cases we can guarantee one directed loop at Vg '■ In the case of a directed 3-cycle G we 
distinguish two cases. Either h)2.i) is fulfilled or in C there exists a 2-path with conditions as in b)2.iii). 
If we have a 2-path which is not a subpath of a directed 3-cycle then we get condition b)2.ii). Hence, $ 
is not bipartite. For the proof of regularity, note, that we consider at each vertex Vg the number of pairs 
of non-adjacent arcs in a realization G. This is the number of all possible arc pairs minus the number of 
adjacent arcs ('^2*^^') — (X)"=i ( 2O ( 2O +X]"=i '^i^i) where J2^=i ( 2O ^^'^ number of all incoming 

arc pairs at each vertex, X]"=i (2') number of all outgoing arc pairs at each vertex and X^iLi '^^^i 

the number of directed 2-paths in a realization G. Hence, the number of non-adjacent arcs is a constant 
value for each realization G. For each of these arc pairs we either set a directed loop or an incoming and 
an outgoing arc at each vertex G V^. For each 2-path in G we set a loop if it is not part of a directed 
3-cycle G = (w^j , , , ) which is an induced subdigraph G = G ({wij , , }). If it is the case it 
exists a realization G' with |GAG'| ~ 6 and GAG' contains exactly 3 different vertices. Hence, we set 
for the 2-path in G with ij < ij> and ij> < ijn with G {1,2,3} the directed arcs {Vg,Vq') and 

{Vg', Vg) and for both other 2-paths in G a directed loop. Generally, we set for all 2-paths in a realization 
an incoming and an outgoing arc at each Vg- The number of 2-paths in each realization is the constant 



value X]r=i ^i^i- Hence, the vertex degree at each vertex is :- 



□ 



In the next section we have to prove that our constructed graphs are strongly connected. This is 
sufficient to prove the reachability of each realization independent of the starting realization. Fig. [5] 
shows an example how the realization G from Fig. [3] can be transformed to the realization G' by a 
sequence of swap operations. 
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G = Go 



swap arcs (5,3) and (2,4) 
with (2,3) and (5,4) 

swap arcs (1,2) and (3,4) 
with (3,2) and (1,4) 

G2 

swap airs (6,5) and (4,7) 
with (6,7) and (4,5) 

0^= G' 



Figure 5: Transforming G from Fig. [3] into G' by a sequence of swap operations. 
3.1 Symmetric differences of two different realizations 

Proposition 3.2. Let S be a graphical sequence and G and G' be two different realizations. If GAG' 
is exactly one weak component and \GAG'\ ^ 6 then there exists in GAG' a vertex-disjoint alternating 
3-walk of type P or Q, where P = (wi, U2, wa, ^4) with (vi,V2), {v3,V4) G A{G) and (u3,W2) G ^(C) and 
Q = {wi,W2,W3, W4) with {wi,W2), (w3, W4) £ ^(C) and (ti'3, 1^2) G A{G). 

Proof. Note that in GAG' an alternating cycle of length two is not possible. Otherwise, there exists an arc 
{u,v) G A{G) n A{G') in contradiction to our assumption that {u,v) S GAG'. The symmetric difference 
GAG' may decompose into a number of alternating cycles {GAG')i. We consider a decomposition into 
the minimum number of such cycles. If one of these alternating cycles {GAG')i contains a vertex- 
disjoint alternating 3-walk P or (5 as claimed, we are done. Otherwise, each vertex is repeated at each 
third step in (GAG')i. Hence, we get the alternating cycles (GAG')i :— (w^ ^ , , , ^ , , , ^ ) 
where (ti^^ , ii^J, (i;^^ , ■^ij), (11^3 , J £ A{G) and (vj^ , J, (u,3 , VjJ, (w^^ , Uig) 6 A{G'). The cycle cannot 
be longer, as the graph induced by GAG'({ui, U2, fa}) is already complete. Since \{GAG')i\ = 6, there 
must be (GAG')j with i ^ j. {GAG')j shares at least one vertex with (GAG')i, because GAG' is 
weakly connected. There must be exactly one = Vj-^^ , since otherwise these two cycles were not arc- 
disjoint. The union of these two cycles is an alternating cycle, in contradiction to the minimality of the 
decomposition. □ 

Note that the above proposition does not assert that the symmetric difference contains P and Q. The 
smallest counter-example are the realizations G = {V,A) and G' = {V,A') with V ~ {wi, U2, wa, W4} and 
A = {(vi,V3), (wa,W2), iv2,V4,), (u4,wi)} and A' = {{vi,V2), (w2,wi), (wa,W4), (w4,wa)}- 

Proposition 3.3. Let S be a graphical sequence and G and G' be two different realizations. If |GAG'| = 
6, then there exist 

a) realizations Go, Gi, G2 with Gq := G, G2 G' and |GiAGi+i| = 4 for i G {0, 1} or 

b) G and G' are different in the orientation of exactly one directed i-cycle. 

Proof. First observe that the symmetric difference is weakly connected whenever |GAG'| = 6. We 
consider the alternating 6-cycle G := GAG' . 
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case 1: C contains at least four different vertices. Assume first that C contains four different vertices. The 
only possibility to realize this scenario is C = {vi,V2, V3,vi,va,v^,vi) with (vi, W2), (1^3, vi), (va, v^) G 
A{G) and (ws, W2), (^4, "^i), (fi, ^^3) G ^(C)- (A permutation of {1,2,3} does not influence the 
result.) We get the alternating vertex-disjoint walk P = (w4, W3, wi, W2)- 

(i) : Assume (v4,U2) ^ A{G). It follows (u4,V2) ^ A{G'). Otherwise, we would get (w4,V2) G GAG' 

in contradiction to our assumption. We set 

Gi := (Go \ {(w4,W3), (wi,W2)}) U {(-y4,W2), (wi,W3)} and 

G2 ;= (Gi \ {{V4,V2), {V3,vi)}) U {(u4,Wl), (W3,W2)}. 

We get G2 = G' and reahzations Go, Gi, G2 with |GiAG-^;^| = 4. 

(ii) : Assume (v4,W2) G A(G). It follows (w4,W2) G A(G'). Otherwise, we would get {v4,V2) G GAG' 

in contradiction to our assumption. We set 

Gi (Go \ {{v4,V2), {v3,vi)}) U {(v4,wi), {v3,V2)} and 

G2 ;= (Gi \ {{V4,V3), {V1,V2)}) U {(W4,W2), {vi,V3)}. 

We get G2 G' and realizations Go, Gi, G2 with |GiAG-^i| = 4. 

We can argue analogously if G contains five our six different vertices. 

case 2: G contains exactly three different vertices. Then G is the alternating cycle G = (wi, U2, t'3, vi, V2-, W3, wi) 
with (?;i,'y2), (^'2,1^3), (■^3,^1) G ^(G) and (^3, W2), (w2, wi), (wi, t'a) G ^(G'). Hence, G and G' are 
different in the orientation of exactly one directed 3-cycle. 

□ 

Lemma 3.4. Let S be a graphical sequence and G and G' be two different realizations. There exist 
realizations Go, Gi, . . . , Gk with Go := G, Gk ■= G' and 

1. |G,AG,+i| =4 or 

2. |G,AG,+i| = 6 

where k < i|GAG'| — 1. In case (2), GjAGi+i consists of a directed 3-cycle and its opposite orientation. 

Proof. We prove the lemma by induction according to the cardinality of the symmetric difference |GAG' | = 
2k. For K := 2 we get |GAG'| = 4. The correctness of our claim follows with Gi :~ G'. For k 3 we get 
a sequence of realizations Gq, Gi, G2 with case a) of Proposition 13.31 In case b) we get a directed 3-cycle 
with its opposite orientation. In both cases it follows k < 2. 

We assume the correctness of our claim for all k < £. Let |GAG'| = 2£-|- 2. We can assume that k > 3. 
Assume further, that the symmetric difference consists of k weakly connected components (GAG')i for 
i G {l,...,fc}. 

Consider first the case that for all these components |(GAG')i| = 6 and that each component contains 
exactly three distinct vertices, then each of them is a directed 3-cycle and its reorientation. We choose 
(GAG')i, perform a 3-cycle reorientation on it, and obtain realization G*. Thus |G*AG'| = 2£ — 4. By 
the induction hypothesis, there are realizations Go = G* ,Gi, . . . ,Gk = G' such that k < i|G*AG'| — 1 < 
5 1 GAG' I — 1. Combining the first 3-cycle reorientation with this sequence of realizations gives the 
desired bound. If there is a component |(GAG')i| = 6 with at least four distinct vertices, we can apply 
Proposition 13.31 case a) to it and handle the remaining components by induction. 

Otherwise, there is a component with |(GAG')j;| > 8. Due to Proposition 13.21 we may assume that 
there is a vertex-disjoint walk P = (wi, W2, W3, W4) with (wi, W2), (^3, U4) G A{G) and (w3,U2) G A[G'). 
Otherwise, there exists Q = (101,^2,^3,^4) with (wi, W2), (103, W4) G A{G') and (1^3,^2) G A{G). In 
that case we can exchange the roles of G and G' and consider G' AG. Clearly, a sequence of realizations 
G' ~ G'q , G^ , . . . , G'j^ = G can be reversed and then fulfills the conditions of the lemma. So from now on 
we work with P. 
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case 1: Assume (wi,U4) e A{G') \A{G). 

This implies (wi, V4) e GAG' . Gi := (G'o\{('yi, W2), (^3, V4)})U{{v3, V2), {vi,V4)} is a realization of 5' 
and it follows IGqAGiI = 4 and |GiAG'| = 2^+2-4 = 2(£-l). Note that after this step, GiAG' may 
consist of several connected components, but each of them has strictly smaller cardinality. There- 
fore, we can apply the induction hypothesis on |Gi AG'|. Thus, we obtain realizations Gi, G2, . . . , Gk 
with Gk ■■= G' and \G^AG^+l\ 4 or |G,,AG,+i| = 6 where fc - 1 < i|GiAG'| - 1. Hence, we get 
the sequence Go, Gi, . . . , Gfe with k = l + i|Gi AG'| - 1 = ^(IGAG'I - 4) < i|GAG'| - 1 which 
fulfills 1. and 2. 

case 2: Assume (wi, W4) e A{G) n A{G'). 

This implies (vi,Vi{) ^ GAG'. Consider an alternating cycle G = {v4,Vi, . . . ,Vj , vi,V2,V3,V4) of 
(GAG')i such that each vertex has in-degree two or out-degree two. Then P is an alternating 
subpath of G with (u^, U4), {vi,Vj) G A{G'). We construct a new alternating cycle G* := (G \ P) U 
{{vijVi)} with length |G*| = |G| — 2. We swap the arcs in G* and get a realization G* of S with 
|GoAG*| = |G*| < 2£ and |G*AG'| = |GAG'| - (|G*| - I) + I < 2£. According to the induction 
hypothesis there exist sequences Gj := G, G\, . . . , G[,^ := G* and Gq :— G* ,G\, . . . , G^^ = G' with 
fci ^ ^|GoAG*| — 1 and fc2 < ^|G*AG'| — 1. We arrange these sequences one after another and get 
a sequence which fulfills 1. and 2. and k = ki+k2 = i|GoAG*| - 1 + ^\G*AG'\ - 1 = ||GAG'| - 1. 

case 3: Assume (wi,U4) £ A{G) \A{G'). 

This implies {vi,V4) G GAG'. The alternating walk P can be extended to an alternating cycle 
G = {vi, V2, W3, W4, W5, . . . , V2t, vi},t>'i using only arcs from GAG' . To construct G, start with P, 
and keep adding alternating arcs until you reach the start vertex vi for the first time. Obviously, 
you will not get stuck before reaching vi. Note that the arc (wi, ^4) does not belong to G. Therefore, 
there exists an alternating sub-cycle G* := GU{(ui,W4)}\P formed by arcs in GAG'. We swap the 
arcs in G* and get a realization G* of S with |GoAG*| = |G*| < 2^ and |G*AG'| = |GAG'|-|G*| < 
21. According to the induction hypothesis there exist sequences Gj :— G, G},...,Gj, := G* 
and Gl := G* ,Gl, . . . ,Gl^ = G' with fci < i|GoAG*| - 1 and fca < i(|GAG'| - |G*|) - 1. 
We arrange these sequences one after another and get a sequence which fulfills 1. and 2. and 
k = ki + k2 = i|GoAG*| - 1 + i(|GAG'| - |G*|) - 1 = i(|GAG'|) - 2. 

case 4: Assume (ui,W4) ^ A{G) U A(G'). This implies (wi,W4) ^ GAG'. It exists the alternating cycle 
G := (P,(ui,U4)) with (ui,W4) ^ ^(G). Gl ;= (Go \ {(ui, ^2), (wa, ^4)}) U {(wg, ^2), ^4)} is a 
reahzation of S and it follows |GoAGi| = 4 and |GiAG'| = 2£ -1- 2 - 2 = 2€. According to the 
induction hypothesis there exist realizations Gi, G2, . . . , G^- with Gk := G' which fulfill 1.) and 2.) 
where fci 1 and ^2 := fc — 1 < i|GiAG'| — 1. Hence, we get the sequence Gq, Gi, . . . ,Gk with 
k = ki + k2 = l + jIGiAG'l - 1 = 5(|GAG'| - 3 1) = i|GAG'| - 1 which fulfills 1. and 2. 

□ 

Corollary 3.5. State graph ^ is a strongly connected directed graph. 

3.2 Random Walks 

A random walk on <t> = (V^, A^) can be described by Algorithm[21 We now require a data structure DS 
containing all pairs of non-adjacent arcs and all directed 2-paths in the current realization. 

Theorem 3.6. Algorithm\^ is a random walk on state graph $ which samples uniformly at random a 
directed graph G' ~ (V, A) as a realization of sequence S for t — >■ 00. 

Proof. Algorithm [5] chooses elements in DS with the same constant probability. For a vertex Vg G 
there exist for all these pairs of arcs in A(G') either incoming and outgoing arcs on Vc in $ or a loop. 
Let d^ := ('^^f "2Y,7=i i"-^)- We get a transition matrix M for $ with ptj = for i,j G A{^),i ^ j, 
Pi] = 1 ^ j)6A(*), i^j} j; for i, J G V^, i j, otherwise we set = 0. Since, $ is a regular, strongly 
connected, symmetrical and non-bipartite directed graph, the distribution of all realizations in a tth step 
converges asymptotically to the uniform distribution. □ 
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Algorithm 2 Sampling realization digraphs 



Input: sequence S, a directed graph G = {V,A) with = (5') Vi G {l,...,n} and E V, a, 

mixing time r. 

Output: A sampled directed graph G' = (V^, A') with {'^^'''^']) = G {1, . . . ,n} and Vi G F. 

1: t := 0, G' G I j initialization 

2: while t < T do 

3: Choose an element p from DS miiformly at random./ /p is a pair of non-adjacent arcs or a directed 
2 -path. 

4; if p is a pair of non-adjacent arcs (vi-^^Vi^), {vi^,Vi^) then 

5: if (ui3,UiJ ^ A{G') then 

6: // Either walk on $ to an adjacent realization G' 

7: Delete {vi^,Vi^),(vi^,Vi^) in A{G'). 

8: Add {v,^,v,^), (wi3,w,;J to A{G'). 

9: else 

10: / / or walk a loop: 'Do nothing ' 

11: end if 

12: else 

13: I Ip is a directed 2-path P ~ (u^j , Vi^ , Vi^) 

14: if ((wj3,Vji) G A{G')) A {{v,^,v^,), {v^,,v^^), iv^,,v,^) ^ ^(G")) A {is > ii) A (^3 > ^2) then 

15: / / Walk on $ to an adjacent realization G' with a reoriented directed S- cycle 

16: Delete {vi^,Vi2), {vi^_,Vi^), {vi^,Vi^) in A{G'). 

17: Add {vi^_,ViJ,{vi^,Vi^_),{vi^,Vi^) to A(G"). 

18: else 

19: //Walk a loop: 'Do nothing' 

20: end if 

21: end if 

22: update data structure DS 

23: t ^t+l 

24: end while 



4 Arc- Swap Sequences 

In this section, we study under which conditions the simple switching algorithm works correctly for 
digraphs. The Markov chain used in the switching algorithm works on the following simpler state graph 
$ = (V^, A^). We define A-^ as follows. 

a) We connect two vertices Vg,Vg' G V-^,G 7^ G' with arcs {Vg,Vg') and {Vg',Vg) if and only if 
I GAG" I = 4 is fulfilled. 

b) We set for each pair of non-adjacent arcs {vij^,Vi^), {vi^,Vi^) G A{G),ij G {1, . . . ,n} a directed loop 
(Vg, Vg) if and only if {v,, , ViJ G A{G) V {vi, , w^J e A{G). 

c) We set one directed loop (Vg, Vg) for all Vg G X^. 

Lemma 4.1. T/ie state digraph <f> = (^^, is non-bipartite, symmetric, and regular. 

Proof. Since each vertex Vg € V^ contains a loop, $ is not bipartite. At each time we set an arc we 
also do this for its opposite direction. Hence, $ is symmetric. The number of incoming and outgoing 
arcs at each Vg equals the number of non-adjacent arcs in G, which is the constant value ('^g'^'") — 

(Er=i ("2 ) + Er=i (2) + ELi a^b?j . Thus, we get the regularity of ¥. □ 
4.1 Characterization of Arc-Swap Sequences 

As shown in Example II. II in the Introduction, <i> decomposes into several components, but we are able to 
characterize sequences S for which strong connectivity is fulfilled in In fact, we will show that there 
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are numerous sequences which only require switching by 2-swaps. In the following we give necessary and 
sufficient conditions allowing to identify such sequences in polynomial running time. 

Definition 4.1. Let S be a graphical sequence and let G = iV^A) be an arbitrary realization. We denote 
a vertex subset V' Q V with \V'\ = 3 as an induced cycle set V' if and only if for each realization 
G* = {V,A*) the induced subdigraph G* {V) is a directed Z-cycle. 

Definition 4.2. Let S be a graphical sequence and G ~ (V', A) an arbitrary realization. We call S an 
arc-swap-sequence if and only if each subset V of vertices with \ V'\ =3 is not an induced cycle set. 

This definition enables us to use a simpler state graph for sampling a realization G for arc-swap- 
sequences. In Theorem 14.51 we will show show that in these cases we have only to switch the ends of two 
non-adjacent arcs. 

Before, we study how to recognize arc-swap sequences efficiently. Clearly, we may not determine all 
realizations to identify a sequence as an arc-swap-sequence. Fortunately, we are able to give a characteri- 
zation of sequences allowing us to identify an arc-swap-sequence in only considering one realized digraph. 
We need a further definition for a special case of symmetric differences. 

Definition 4.3. Let S be a graphical sequence and G = (V, A) and G* ~ {V,A*) arbitrary realizations. 
We call GAG* simple symmetric cycle if and only if each vertex v G y(GAG*) possesses vertex in-degree 
'^GAG' (^) — ^ ^^'^ vertex out-degree d'^^^^^-^^, < 2, and if GAG* is an alternating directed cycle. 

Note that the alternating directed cycle Gi in Fig. U] is not a simple symmetric cycle, because 
(ij^(4) = 4. Cycle Gi decomposes into two simple symmetric cycles G[ = {wi, W2, "^3, "^4, wi} and 

G" = {V2,V3,V5,V4,V2}- 

Theorem 4.2. A graphical sequence S is an arc-swap-sequence if and only if for any realization G = 
{y, A) the following property is true: 

For each induced, directed 3-cycle G (V) of G there exists a realization G* = (V,A*) so that GAG* is a 
simple symmetric cycle and that the induced subdigraph G* {V) is not a directed S-cycle. 

Proof. =>: Let S* be a graphical arc-swap sequence and G = {V,A) be an arbitrary realization. With 
Definition 14.21 it follows that each subset V' C V with \V'\ = 3 is not an induced cycle set. Hence, there 
exists for each induced, directed 3-cycle G {V') of G a realization G' = {V, A') with symmetric difference 
GAG' where the induced subdigraph G' (V) is not a directed cycle. If the symmetric difference GAG' is 
not a simple symmetric cycle we delete as long alternating cycles in GAG' as we get an directed alternating 
cycle G* where each vertex in G* has at most vertex in-degree two and at most vertex out-degree two. 
Furthermore, G* shall contain at least one arc {v,v') ^ V' x V' . This is possible, because GAG' contains 
at least one such arc. On the other hand the alternating cycle G* does not contain all possible six of such 
arcs. Otherwise, the induced subdigraph G' (V) is a directed cycle. Now, we construct the realization 
G* = {V, A*) with A* := (A(G) \ {A{G*) n A{G))) U {A{G*) n A{G')). It follows GAG* = G* is a simple 
symmetric difference. 

<;=: Let G be any realization of sequence S. We only have to consider 3-tuples of vertices V inducing 
directed 3-cycles in G. With our assumption there exists for each V' a realization G* so that G* {V') 
is not a directed 3-cycle. Hence, we find for each subset V' C V oi vertices with = 3 a realization 
G* — {V, A), so that the induced subdigraph G* (V) is not a directed 3-cycle. We conclude that S is an 
arc-swap sequence. □ 

This characterization allows us to give a simple polynomial-time algorithm to recognize arc-swap- 
sequences. All we have to do is to check for each induced 3-cycle of the given realization, if it forms 
an induced cycle set. Therefore, we check for each arc (w, w) in an induced 3-cycle whether there is an 
alternating walk from v to w (not using arc {v,w)) which does not include all five remaining arcs of the 
3-cycle and its reorientation. Moreover, each node on this walk has at most in-degree 2 and at most 
out-degree 2. Such an alternating walk can be found in linear time by using a reduction to an /-factor 
problem in a bipartite graph. In this graph we search for an undirected alternating path by growing 
alternating trees (similar to matching algorithms in bipartite graphs, no complications with blossoms will 
occur), see for example |Sch03| . The trick to ensure that not all five arcs will appear in the alternating 
cycle is to iterate over these five arcs and exclude exactly one of them from the alternating path search 
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between v and w. Of course, this loop stops as soon as one alternating path is found. Otherwise, no such 
alternating path exists. As mentioned in the Introduction, a linear-time recognition is possible with a 
parallel Havel-Hakimi algorithm of LaMar |LaM09| . 

Next, we are going to prove that <I> is strongly connected for arc-swap-sequences. The structure of 
the proof is similar to the case of $, but technically slightly more involved. 

Lemma 4.3. Let S be a graphical arc- swap- sequence and G and G* be two different realizations. Assume 
that V' := {vi,V2,V3} C V such that G{V') is an induced directed 3-cycle but G*{V') is not an induced 
directed 3-cycle. Moreover, assume that GAG* is a simple symmetric cycle. Then there are realizations 
G'o,Gi,...,Gfc with Go ■.= G,Gk := G* , \G,AG^+i\ ^ A and k < ^\GAG*\. 

Proof. We prove this lemma by induction on the cardinality of GAG*. The base case \GAG*\ = 4 is 
trivial. Consider next the case \GAG* \ = 6. We distinguish between two subcases. 

case a) GAG* consists of at least four different vertices. 

By Proposition 13.31 case a), there are realizations G ~ Gq,Gi,G2 = G* with \GiAGi+i \ = 4. 

case b) GAG* consists of exactly three vertices V4,V5,ve. 

Observe that GAG* contains at least one arc from G{V') or its reorientation but not all three ver- 
tices V as otherwise G*{V') would be an induced 3-cycle. In fact, it turns out that GAG* contains 
exactly one arc, say (w2, ws), from G(V') and its opposite arc (wa, W2), because GAG* is the directed 
alternating cycle C := (i'2, W3, W4, ^2, W3, W4, W2) with V4 ^ vi. We have two subcases. Assume first 
that {vi,V4) ^ A{G) n A{G*). So we can swap the directed alternating cycle (wi, W2, ^3, W4, vi) in a 
single step. We then obtain the directed alternating 6-cycle (y2,V3,Vi,V2,vi.,V4,V2) which consists 
of four different vertices. By Proposition 13.31 case a) , we can swap the arcs of this cycle in two 
steps, thus in total in three steps as claimed. Otherwise, {vi,V4) G A{G) fl A(G*). Then we obtain 
the directed alternating cycle (wi, W4, W2, W3, wi) which can be swapped in a single step. By that, 
we obtain a new cycle (ui, V3, W4, U2, W3, W4, wi) which consists of four different vertices. By Propo- 
sition [2]31 case a), we can swap the arcs of this cycle in two steps, thus in total in three steps as 
claimed. 

For the induction step, let us consider |GAG*| = 2£ -\- 2 > 8. Then GAG* contains between one 
and five arcs from G{V') and its reorientation. By Proposition 13.21 there is a vertex-disjoint alternating 
directed walk P = {wi,w2,w3,w4) in GAG* with {wi,w2) G A{G)\A{G*) or {wi,w2) G A{G*)\A{G). 

Suppose that P contains no arc from G{V') and its reorientation. We consider the case {wi,W2) S 
A{G) \ A{G*). If {wi,w4) G A{G) n A(G*), then we consider G = (GAG*) U {{wi,W4)} \ P. We 
swap the arcs of G and obtain as realization G**. Clearly, GAG** contains an arc from G{V') or 
its reorientation, and is a simple symmetric cycle. As |GAG**| = 2£, we can apply the induction 
hypothesis. We obtain a sequence of realizations G = Gq,Gi, . . . ,Gk = G** with \GiAGi+i\ = 4 and 
k < i I GAG** I < i (I GAG* I - 2). Finally, we apply a last swap on the cycle {wi,W2,W3,W4,wi) and 
thereby transform G** to G*. In total, the number of swap operations is A: < i|GAG*|. 

The case {wi,W4) ^ A{G) D A{G*) is similar. This time, we start with a single swap on the cy- 
cle (wi,W2,W3,W4,wi) and afterwards apply induction to the remaining cycle. We can treat the case 
{wi,W2) G A{G*) \ A{G) analogously. Thus we can exclude the existence of any vertex-disjoint directed 
alternating 3-walk which does not contain at least one arc from G{V') and its reorientation. 

It remains to consider the case that there is a vertex-disjoint directed alternating 3-walk P = 
{wi,W2,W3,W4) in GAG* with (^1,^2) G A{G) \A{G*) or {■wi,W2) G A{G*) \A{G) but at least one arc 
of P is from G{V') and its reorientation, say {vi, V2). 

Recall that GAG* contains between one and five arcs from G{V') and its reorientation. We distinguish 
between three cases: 

case I: GAG* contains exactly one of these arcs, say {vi,V2) G A{G) \ A{G*). (The case that {v2,vi) G 
A{G*) \ A{G) is this special arc can be treated analogously.) 

We claim that the cycle GAG* must have the form (ui, 1)2, W4, W5, we, W4, W5, wg, wi). Note that 
V4,V5,vg are repeated every third step, as otherwise we would obtain a vertex-disjoint alternat- 
ing cycle as excluded above. The cycle cannot be longer than eight, since then we would either 
obtain a vertex-disjoint 3-walk (z;4, ^5, wg, W7), also excluded above, or if V4 = vj we would violate 
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simplicity of the symmetric difference, ft might be that = v^, but Vi,ve ^ as otherwise 
the symmetric difference would contain more than one arc from G{V') and its reorientation, ff 
{vi,V4) S A{G) n A{G*) there is the alternating directed 4-cycle {vi,V4,V5,ve,vi) which can be 
swapped. In the remaining 6-cycle the arc (wi,W2) is contained, so the induction hypothesis can 
be applied. Otherwise, if {vi,V4) ^ A{G) D A{G*), we first apply the induction hypothesis to the 
6-cycle (ui, U2, i'4, i^s, we, W4, ui), and afterwards we swap the remaining 4-cycle (ui, V4,V5,ve,vi). 

case ff: GAG* contains exactly two of these arcs. 

Suppose first that these two arcs are adjacent, say (vi, W2), (^3, V2). Consider the following arcs 
iv3-,V4)^ {v^yVi), (v5,ve) along the symmetric difference. Now = V2 as otherwise there is an 
alternating directed walk (i'2, "^3, f4, 1's)- Depending whether (y5,V2) € A{G) Ci A{G*) or not, we 
can either swap the alternating 4-cycle (^3, W4, W5, W2, ^3) or the remaining part of the symmetric 
difference together with (w.g, V2) by the induction hypothesis. Moreover, vq = W3, as otherwise there 
would be the vertex-disjoint alternating directed 3-walk (v3,V4,vr,,ve) excluded above. But then 
{v3, V2) is also in the symmetric difference, a contradiction. Thus, the two arcs from G{V') and its 
reorientation are not adjacent. Then, there are at least two other arcs between them (otherwise 
the one arc between them would also be from G{V') and its reorientation). By our assumption, 
there is a vertex-disjoint alternating directed 3-walk P = (u>i, u'2, 1^3, ^4) with at least one arc from 
G{V') and its reorientation. In our scenario it must be exactly one such arc. Depending whether 
{wi,W4) G A(G) n A(G*) or not, we can either swap the alternating 4-cycle {wi,W2,W3,W4,wi) or 
the remaining part of the symmetric difference together with {wi,W4) by the induction hypothesis. 

case Iff: GAG* contains between three and five of these arcs. 

Suppose first all of them follow consecutively on the alternating directed cycle. Consider the last 
two of these arcs, and append the next arc which must end in a vertex V4 ^ V' . Then we have a 
vertex-disjoint alternating directed 3-walk which contains two arcs from G{V') and its reorientation, 
and the remaining part of the symmetric difference has also such an arc. Thus we can apply the 
induction hypothesis and are done. Otherwise the three to five arcs from G{V') and its reorientation 
are separated. So no alternating directed 3-walk may contain all of them, in particular not P. We 
can proceed as in case ff). 

□ 

Proposition 4.4. Let S be a graphical arc-swap-sequence and G and G' be two different realizations. If 
\GAG'\ = 6 and GAG' consists of exactly three vertices V' :~ {i'i,W2,W3}, then there exist realizations 
Go,Gi,...,Gfc with Go ■.^G,Gk := G', \G^AGi+l \ =4 and k <2n + 2. 

Proof. Since S is an arc-swap-sequence. Theorem 14.21 implies the existence of a realization G* such that 
GAG* is a simple symmetric cycle and G* (V) is not a directed 3-cycle. By Lemma 14.31 there are 
realizations Go,Gi, . . . ,Gk' := G* with |GiAGi+i| ~ 4 and k' < ^\GAG*\ < n, since GAG* is simple. 
Moreover, we have |G*AG'| < |GAG*| + 4 < 27i + 4 since G and G' differ only in their orientation of the 
3-cycle induced by V . The symmetric difference G* AG' is not necessarily a simple symmetric cycle, but 
can be decomposed into simple symmetric cycles, each containing at least one arc from G{V') and its 
reorientation. On each of these simple symmetric cycles we apply our auxiliary Lemma 14.31 We obtain 
a sequence G* := Gg, . . . , G'j.,, := G' with \GiAGi+i \ = 4 and k" < n + 2. Combining both sequences we 
obtain a sequence with k = k' + k" < 2n + 2. □ 

Lemma 4.5. Let S be a graphical arc-swap-sequence, and G and G' be two different realizations. 
Then there exist realizations Gq, Gi, . . . , Gfe with Gq '.— G, Gk '■— G' and \GiAGi+i\ = 4, where 
k < (i|GAG'| - f) • (71 + 1). 

Proof. We prove the lemma by induction according to the cardinality of the symmetric difference |GAG' | = 
2k. For K := 2 we get |GAG'| — 4. The correctness of our claim follows with Gi := G'. 

For K := 3 we distinguish two cases, ff GAG' consists of exactly three vertices, then by Proposition 14. 41 
we get a sequence of realizations Gq, Gi, . . . ,Gk and k < 2n + 2 — 2{n -f 1), as claimed. Otherwise, the 
symmetric difference GAG' consists of more than three vertices. By Proposition 13.31 case a), there are 
realizations Gq, Gi, G2 = G'. 

We assume the correctness of our induction hypothesis for all k < £. Let |GAG'| = 2£ + 2. We can 
assume that k > 3. Suppose first that the symmetric difference GAG' decomposes into t simple symmetric 
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cycles \{GAG')i\ = 6. Suppose further that all these (GAG"); consist of exactly three vertices. Clearly, 
|GAG'| = 6t. We apply our Proposition 14. 41 to each of these t cycles one after another and get a sequence 
of realizations Go,Gi,...,Gfe = G' with k < 2t{n + 1) < (3i- l)(n+ 1) = (i|GAG'| - l) ■ (n + 1). 

Otherwise, there is a (GAG')i which contains at least four vertices. Swapping the arcs in (GAG')i 
leads to a realization G*. By Proposition 13.31 there are realizations G = Go,Gi,G2 = G* with 
|GiAGi+i| ~ 4. We can apply the induction hypothesis on the remaining part of the symmetric dif- 
ference. Obviously, we obtain the desired bound in this case. 

It remains the case that there exists a simple symmetric cycle (GAG')i of GAG' with |(GAG')i| 6. 
If |(GAG')i| = 4, we use a single swap on (GAG')i and obtain a realization G* , where |G*AG'| = 
I GAG' I — 4. By the induction hypothesis, there is a sequence of realizations G* = Gi,G2 ■ ■ ■ ,Gk = 
G' with fc - 1 < (i|G*AG'| - 1) • (n + 1) = (i|GAG'| -3) • (n + 1). Otherwise, |(GAG')i| > 8. 
Using Proposition 13.21 we may assume that there exists a vertex-disjoint directed alternating walk P = 
(wi, U2, W3, W4) in (GAG')i with {vi,V2), (w3,W4) <S A{G) and (w3,W2) & A{G'), for the same reasons as in 
the proof of Lemma 13.41 

case 1: Assume (wi,U4) 6 A{G') \A{G). 

This implies {vi, V4) G GAG' . Gi := (G \ {{vi,V2), (fa, V4)}) U {(fa, i'2), (^"1, W4)} is a realization of 
S and it follows |GAGi| = 4 and |GiAG'| = 21 + 2 - 4 = 2{i - 1). Therefore, we can apply the 
induction hypothesis on |GiAG'|. Thus, we obtain realizations Gi, G2, . . . ,Gk with Gk G' and 
|G,AG,+i| =4. where fc - 1 < (i|GiAG'| - l) • (n + 1) = (i|GAG'| -3) • (n + l). 

case 2: Assume (wi,U4) £ A{G) nA{G'). 

This implies (ui,U4) ^ GAG'. We construct a new alternating cycle C* :— ((GAG')i\P)U{(?ji, U4)} 
with length |G*| = |(GAG')i| — 2. We swap the arcs in G* and get a realization G* of S with 
|GoAG*| = |G*| < 2£ and |G*AG'| = |GAG'| - (|G*| - 1) + 1 < 2£. According to the induction 
hypothesis there exist sequences Gq := G, Gj, . . . , G^^ := G* and Gq := G*, G^, . . . , G^^ = G' with 
h < (i|GoAG*| - 1) • (n + 1) and k2 < (i|G*AG'| - l) • (n + 1). We arrange these sequences one 
after another and get a sequence with = fci + fc2 = (i|GoAG*| - 1 + ^|G*AG'| - l) • (71 + 1) = 
(i|GAG'|-l)-(n + l). 

case 3: Assume (wi, U4) 6 A{G) \ A{G'). 

This implies (ui,W4) G GAG'. Note that the arc {vi,V4) does not belong to (GAG')i. There- 
fore, there exists an alternating sub-cycle G* := (GAG')i U {{vi,V4)} \ P formed by arcs in 
GAG'. We swap the arcs in G* and get a reahzation G* of S with |GoAG*| = |G*| < 2£ and 
G*AG'| = I GAG' I — |G*| < 21. According to the induction hypothesis there exist sequences 
Gl := G,G\,...,Gl^ := G* and G^ := G*, G?, . . . , G^.^ = G' with fci < (i|GoAG*| - l) ■ {n + l) 
and k2 < (i(|GAG'| — |G*|) — l) ■ (n-|- 1). We arrange these sequences one after another and get a 
sequence with k = fci-hfe = (i|GoAG*| - 1 + i(|GAG'| - |G*|) - l) -(n-hl) = (i(|GAG'|) - 2) • 
(n + l). 

case 4: Assume {vi,V4) ^ A{G) U A{G'). This implies {vi,V4) ^ GAG'. It exists the alternating cycle 
G := (P,(ui,U4)) with {vi,V4) i A{G). Gi := (Go \ {(wi, W2), ("a, W4)}) U {(wa, W2), (wi, W4)} is 
a reahzation of S and it follows |GoAGi| = 4 and |GiAG'| = 2^ -I- 2 - 2 = 2£. According to 
the induction hypothesis there exist realizations Gi,G2,...,Gfc with Gfe G' where k\ := 1 
and k2 := fc — 1 < (i|GiAG'| — l) • (n + 1). Hence, we get the sequence Go,Gi,...,Gfe with 
fc = fci + fc2 = 1 + (^IGiAG'l - 1) • (n + 1) < (i|GAG'| - l) • (n + 1). 

□ 

Corollary 4.6. State graph $ is a strongly connected directed graph if and only if a given sequence S is 
an arc-swap-sequence. 

An arc-swap-sequence implies the connectedness of the simple realization graph Therefore, for 
such sequences we are able to make random walks on the simple state graph <I> which can be implemented 
easily. We simplify the random walk Algorithm [5] for arc-swap-sequences in using realization graph $. 
Hence, our data structure DS only contains pairs of non-adjacent arcs. We can ignore lines 12 to 20 in 
Algorithmic] We denote this modified algorithm as the Arc-Swap-Realization-Sample Algorithm 3. 
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Theorem 4.7. Algorithm 3 is a random walk on the state graph <I> which uniformly samples a directed 
graph G' = (V, A) as a realization of an arc-swap-sequence S for r oo. 

Proof. Algorithm 3 chooses all elements in DS with the same constant probability. For a vertex Vq G 
there exist for all these pairs of arcs in A{G') either incoming and outgoing arcs on Vc G Vj- or a 
loop. We get a transition matrix AI for <i> with pij = ^ for i,j S A-^, i ^ j , pij = 1 — } 1 ^'^^ 

i, j G V^, i = j, otherwise we set pij = where d := (''^2'^") —2 X^ILi ("2') ~^7=i Since, $ is a regular, 
strongly connected, symmetric, and non-bipartite directed graph, the distribution of all realizations in a 
tth step converges asymptotically to the uniform distribution, see Lovasz |Lov96| . □ 

4.2 Practical Insights And Applications 

As mentioned in the Introduction, many "practitioners" use the switching algorithm for the purpose of 
network analysis, regardless whether the corresponding degree sequence is an arc-swap-sequence or not. 
In this section we would like to discuss under which circumstances this common practice can be well 
justified and when it may lead to wrong conclusions. 

What would happen if we sample using the state graph $ for a sequence S which is not an arc- 
swap-sequence? Clearly, we get the insight that $ has several connected components, but as we will see 
$ consists of at most 2^~^ isomorphic components containing exactly the same realizations up to the 
orientation of directed 3-cycles each consisting of an induced cycle set V' . Fortunately, we can identify 
all induced cycle sets using our results in Theorem 14.21 bv only considering an arbitrary realization G. 

Proposition 4.8. Let S be a graphical sequence which is not an arc- swap- sequence and has at least two 
different induced cycle sets V' and V" . Then it follows V' n V" ~ 0. 

Proof. Without loss of generality we can label the vertices in V' with v[,v'2,v'^ and in V" with v": ^2 , W3 . 
Let G be a realization where {(vi, U2), iv'2,v'^), (wg, v[), (v", ), ("2 , "3 ), (wg , v'{)} C A{G). We distinguish 
between two cases. 

a) : Assume \V' V"\ = 1 where v[ = v'{. If it exists arc (wg'jUg) S ^(G) we find the alternating 4- 

cycle (wg , Ug, vj^, ^2 1 which implies a new realization G* where G* {V) is not a directed cycle in 
contradiction to our assumption that V' is an induced cycle set. Hence, it follows (wg,Wg) ^ A{G). 
In this case we find the alternating cycle {v'^ v'l, v'2, v'^, v'^) . 

b) : Assume \V' n V"\ = \{v[,v'2}\ = 2 where v'l = v'{ and v'2 = v!^. If arc (ug , Wg) ^ A{G) exists we find 

the alternating 4-cycle {v[,v'^,v'^,V2,v[) which implies a new realization G* where G* (V) is not a 
directed cycle in contradiction to our assumption that V' is an induced cycle set. Hence, it follows 
(wg , I'g) G A{G). In this case we find the alternating cycle (ug , Vg, v[, V2, Ug ). 

□ 

As the induced 3-cycles which appear in every realization are vertex-disjoint, we can reduce the in- 
and out-degrees of all vertices in these cycles by one, and obtain a new degree sequence which must be 
an arc-swap sequence. 

Theorem 4.9. Let S be a sequence. Then the state graph $ consists of at most 2L 3 J isomorphic 
components. 

Proof. We assume that S is not an arc-swap-sequence, otherwise we apply Theorem l4.6l and get a strongly 
connected digraph $. With Proposition 14.81 it follows the existence of at most [-^J induced cycle sets 
for S. Consider all realizations G^ possessing a fixed orientation of these induced 3-cycles which implies 
G-' (Vi) = G^ {Vi) for all such realizations. We pick out one of these orientation scenarios and consider 
the symmetric difference G-' AG^ of two such realizations. Since, all induced 3-cycles are identical in 
G^ and G^ , we can delete each arc of these induced cycle sets V' and get the reduced graphs G^ and 
Gl . Both are realizations of an arc-swap-sequence S". Applying Theorem 14.51 we obtain, that there exist 
realizations Go := G^, . . . ,Gk G{ |GiAGi+i| = 4 and k < \GlAGl \. Hence, each induced subdigraph 
$ (^{Vqj IVqj G V-^ and G^ is a realization for one fixed orientation scenario}^ is strongly connected. On 
the other hand, we get for each fixed orientation scenario exactly the same realizations G-'. Since, all 
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induced 3-cycles are isomorphic, it follows that all realizations which are only different in the orientation of 
such directed 3-cycles are isomorphic. By Theorem l4.2l there does not exist an alternating cycle destroying 
an induced 3-cycle. Hence, the state graph <I> consists of exactly 2*^ strongly connected isomorphic 
components where k is the number of induced cycle sets V' . □ 

Applications in Network Analysis Since the switching algorithm samples only in one single compo- 
nent of (f>, one has to be careful to get the correct estimations for certain network statistics. For network 
statistics on unlabeled graphs, it suffices to sample in a single component which reduces the size of by 
a factor 2'"', the number of components in $, where k is the number of induced cycle sets of the prescribed 
degree sequence. Examples where this approach is feasible are network statistics like the average diameter 
or the motif content over all realizations. 

For labelled graphs, however, the random walk on systematically over- and under-samples the 
probability that an arc is present. Suppose that the random walk starts with a realization G = {V,A). 
If an arc (wi,W2) £ ^(G) belongs to an induced cycle set, it appears with probability 1 in all realizations 
of the random walk. The opposite arc {v2tVi) ^ ^(G), will never occur. In an unbiased sampling over 
all realizations, each of these arcs, however, occurs with probability 1/2. All other arcs occur with the 
same probability in a single component of V-^ as in the whole state graph. This observation can be used 
to compute correct probabilities for all arcs. 

5 Concluding Remarks 

In this paper, we have presented Markov chains for sampling uniformly at random undirected and directed 
graphs with a prescribed degree sequence. The key open problem remains to analyze whether these 
Markov chains are rapidly mixing or not. 
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Appendix 



A Further Examples Where Switching Fails 

As we have seen in Example 11.11 the switching algorithm will fail in general. Here we give further non- 
trivial classes of graphs where it also fails. All problematic instances are realizations which are different 
in at least one directed 3-cycle but not all of them are not changeable with alternating 4-cycles. 

Consider the following Figures [5] and [T] Both examples cannot be changed to a realization which is 
only different in the orientation of the directed 3-cycle by a sequence of alternating 4-cycles. 



ncident to each vertex 




arbitrary digraph 



Figure 6: All vertices in a 3-cycle are incident in one direction with vertices in an arbitrary subdigraph. 




Figure 7: All vertices in a 3-cycle are incident in both directions with a directed clique. An independent 
set of vertices is arbitrarily incident with the directed clique. 
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